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Abstract: We uncover an unexpected connection between the physics of loop integrals 
and the mathematics of spline functions. One loop integrands are Laplace transforms of 
splines. This clarifies the geometry of the associated loop integrals, since a n-node spline 
*^ has support on a n-vertex polyhedral cone. One-loop integrals are integrals of splines 

on a hyperbolic slice of the cone, yielding polytopes in AdS space. Splines thus give a 
geometrical counterpart to the rational function identities at the level of the integrand. 
Spline technology also allows for a clear, simple, algebraic decomposition of higher point 
loop integrals in lower dimensional kinematics in terms of lower point integrals - e.g. an 
hexagon integral in 2d kinematics can be written as a sum of scalar boxes. Higher loops 
can also be understood directly in terms of splines - they map onto spline convolutions, 
leading to an intriguing representation in terms of hyperbolic simplices integrated over other 
hyperbolic simplices. We finish with speculations on the interpretation of one-loop integrals 
as partition functions, inspired by the use of splines in counting points in polytopes. 
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1 Introduction 

Scattering amplitudes are beautiful multi-faceted objects. Nowhere does this show more 
clearly than in the maximally supersymmetric SU(N) Yang-Mills theory [1], especially in 
the planar (N — > oo) sector [2, 3]. There, new perspectives on scattering amplitudes as 
expectation values of null-polygon shaped Wilson loops [4] eventually led to the discovery 
of dual conformal symmetry [5-7] and the consequent infinite symmetry algebra known 
as the Yangian [8]. Relatedly, scattering amplitudes have been shown to be intimately 
connected to correlation functions of stress-tensor multiplets [9, 10]. Concurrently, another 
intriguing picture seems to be emerging connecting the Yangian invariants to the algebra 
and combinatorics of the positive part of the Grassmannian - the space of k x n-dimensional 
planes which satisfy certain positivity conditions [11-15]. 
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The wealth of riches does not stop here. Mellin space [16-18] is the natural place to 
discuss conformal correlation functions, and has led not only to major technical progress 
but also to new perspectives on bulk locality and the flat-space limit of AdS scattering 
amplitudes [19-25]. Whether correlation functions are computed from an AdS dual [26, 27] 
or directly in flat-space [28] , a perturbative formulation ensures that their Mellin amplitudes 
take a simple form highly reminiscent of scattering amplitudes themselves. Mellin space 
also sheds light on the various identities that have been uncovered relating different kinds of 
conformal loop integrals. For instance, in [28] a simple formula was written down connecting 
the scalar double box integral to the d = 6 hexagon. Many other such identities can be 
written down, owing to the particular fact that Mellin amplitudes of flat space correlation 
functions seem to satisfy Feynman rules. In particular, Mellin amplitudes for position space 
tree-level diagrams take a factorized form, which implies that the corresponding position 
space formulae take the form of convolutions of contact interactions. The above mentioned 
identity is the simplest case for this, with a 3 — > 3 exchange diagram being written in 
terms of a 6-pt contact interaction 1 . Analogously, one could write a ladder diagram with 
n-boxes in terms of a simple (n — l)-fold integral of a 4 + 2n contact interaction. It is 
therefore especially important to understand the latter integrals. They can be thought of 
as 4 + 2n-gon integrals in as many dimensions, dimensionally reduced to 4 dimensions, and 
a systematic understanding of these was one of the main motivations for this paper. 

Curiously, the answer to this problem arose from an examination of yet another per- 
spective on loop integrals, which connects scattering amplitudes to convex geometry. This 
has so far been seen to some extent from two different sides, both at level of the integrand 
[29-32], as well as the integral [33]. In the first, the BCFW [34] recursion relation acquires 
a remarkable interpretation, that of a simplicial decomposition of volumes of polytopes in 
projective space. In the second, the integrals were seen to compute volumes of polytopes, 
or boundaries thereof, in hyperbolic space. 

It is with this last approach that this note is concerned with. We examine loop inte- 
grals, conformal and otherwise, from a new perspective, connecting loop integrals to the 
mathematics of splines 2 - a special class of functions better known for its uses in approxi- 
mation theory! The splines have positivity and convexity built in and provide the missing 
link in connecting the integrand and the integral. There are several kinds of splines. The 
ones we will concern ourselves with are sometimes called multivariate truncated powers, 
and they are piecewise polynomials with support on polyhedral cones, whose vertices cor- 
respond to vectors in projective space. The loop integral integrates over the spline with a 
gaussian measure, reducing the integral from the cone to a hyperbolic slice thereof, leading 
to the appearance of polyhedra in AdS. It is important to emphasize that our methods 
are in no way restricted to N = 4 SYM - they are universal and apply to any quantum 
field theory, conformal or otherwise, and are really statements on integrals of products of 
propagators. In particular, this holds even for the calculation of Witten diagrams in AdS 
space! In this way we will find a beautiful generalization of the star-triangle relation, where 

1 In the dual momentum space, these are scalar double box and hexagon integrals, respectively 
2 The literature on splines is vast and we cannot do justice to it here. We have found the book [35] and 
articles [36, 37] especially useful and a large list of references can be found there. 
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many-legged AdS stars get glued up into AdS polyhedra. Crucial to our approach will be 
the use of the embedding space formalism: the spline naturally lives in this space and it is 
from there that the polyhedra descend. 

Let us discuss the outline of this paper. After some preliminaries in section 2 we 
move on to section 3, where we give a new point of view on the old familiar scalar box 
integral. This allows us to introduce the spline in a simple case which captures most of the 
important geometric features of more complicated integrals. Section 4 greatly generalizes 
these results, showing that hyperbolic polyhedra hide in disguise in a much broader class 
of integrals, among which we include contact interactions in AdS. Accordingly we get a 
nice map from AdS stars into AdS polyhedra. In section 5 we emphasize that the spline 
is nothing but the Laplace transform of the one-loop integrand. This allows one to trade 
identities of rational functions for geometric ones, and vice versa. We examine in particular 
how the boundary of a four-simplex can be alternatively be given as a pentagon integral or 
as a sum of boxes, clarifying the simple representations found in [32]. The following section 
deals with the calculation of splines which have large numbers of nodes. What this means 
is that geometrically these correspond to splines associated with non-simplicial cones; from 
the point of view of the integrand it corresponds to considering loop integrals involving a 
number of external points greater than d + 2. We present two methods for computing the 
spline in these cases, one intuitively obvious in terms of partial fraction decomposition of 
the integrand; and another, more powerful in terms of a residue computation. As a simple 
application, we show how to predict the result for a box integral in Id kinematics from 
a purely geometric perspective. We comment on more advanced applications, but leave 
detailed calculations for a future publication [38]. Next we move on to considerations of 
higher loop integrals. We begin by showing that one-loop integrals can be written in several 
different equivalent ways, and how one spline can actually be decomposed into two or more. 
Then we perform an explicit two-loop computation which shows a very similar structure. 
In both cases we show that the spline picture takes the form of integrals of hyperbolic 
simplices over other simplices. This shows that the spline picture is not intrinsic to one- 
loop computations, although the physical meaning of this picture is far from obvious. We 
finish with a discussion and some speculations, in particular on the idea of the spline as 
a charge degeneracy. The paper is complemented by two appendices with some technical 
details on the two-loop calculation and the link between the spline we use and Carlson's 
Dirichlet B-splines [36]. 

2 Set-up 

2.1 Ambient space 

In the calculations to follow, it is crucial that they should be seen in the light of the 
embedding or ambient space formalism [39-42] . In this approach we think of d-dimensional 
Euclidean space as a slice of the projective light-cone in d + 2 dimensions. The SO(d + 
1,1) invariance group acts linearly on d + 2 dimensional vectors and induces a non-linear 
representation of the conformal group in d dimensions. We will consider null vectors in 
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d + 2 dimensions which we will invariably denote by capital letters P and Q, the latter 
being used for dummy variables. These are projective null vectors which means we have 



P M P M = -P+P- + P^P^ = 0, P ~ XP (2.1) 

where M, n are d + 2 and dimensional indices respectively. 

We define a reference vector /, or infinity vector, which explicitly breaks conformal 
SO(d+ 1, 1) invariance. The vector P M /{—P ■ I) parameterizes d-dimensional flat space. 
We choose the gauge — P ■ I = ^?-P + , such that 



2 

pM 



2(l,x 2 ,xn, (2.2) 



-PI 

By contracting two independent P vectors we obtain 

P 'i - (-p'^-p] ■ I) = (x ' ~ x ' )2 - < 2 ' 3 ' 

The factors Pi ■ I drop out of conformally invariant expressions, and we can in any case 
recover them by demanding homogeneity zero in rescalings of the Pi. 

So much for null vectors. We will also need to consider massive vectors, satisfying 
X 2 = — 1 . These parameterize a d + 1-dimensional AdS space and can be written 

X = — (l,x 2 + x 2 ,x fl ) . (2.4) 

Xq 

The coordinates xo and x^ are simply those in the usual parameterization of the Poincare 
patch of AdS, 

dsLs = M±J%L±. ( 2 .5) 
x 

2.2 Dual conformal symmetry 

As mentioned in the introduction, N = 4 SYM has a dual conformal symmetry. We will 
not focus particularly on N = 4 in the following, but we will be using the idea that the 
same integral can be thought of in two dual manners. What this means in practice is that 
if one does the replacements 

Pi -)■ Xi - xj (2.6) 

in a momentum space planar loop integral, the result looks like the calculation of a position 
space correlation function. This is illustrated for the double box integral in figure 1. 
In practice we are replacing propagators as 

Pi ~^ (xi - x i+i ) 2 Pij ^ 2 ' 7 ^ 

It is from this position space perspective that we will be working. We are not following 
the convention that (xi — Xj+i) 2 = 0. Some of the integrals we will consider will have no 
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Figure 1. A fully massive two- loop scalar box integral (black) and its position space counterpart 
(blue). 

interpretation from the momentum space point of view, since we will consider propagator 
factors such as l/(xi — X2) 2 ^ with A ^ 1. These do appear however in AdS/C¥T cal- 
culations. We will keep an attentive eye regarding applications to momentum space loop 
integrals. 

We shall also consider massive momentum space loop integrals. To what does this 
correspond in the dual space? Well, such a propagator is given by 

1 1 1 (2-8) 



p 2 + m 2 (x — y) 2 + m 2 m (x — y) 2 + m 2 
This suggests we identify m as the xq coordinate of an AdS vector X, since we have e.g. 

1 x 



2 ' 



(2.9) 



The mass then becomes the coordinate xo in AdS. In this note we will happily consider 
massless vectors P or massive vectors X as our external data. When a dual momentum 
space perspective is warranted, these correspond to loop integrals with massive or massless 
propagators. In practice, throughout this paper we will be thinking of the integrals as 
position space tree-level correlators. Nevertheless, we will keep the convention of denoting 
the number of integrals involved by the loop-order - which equals the number of interaction 
vertices in the position space picture. 

3 Warm-up: the box and the spline 

As a warm-up to greater things to come, let us consider the familiar scalar box integral. 
In the dual conformal picture, this is the same as a position space contact interaction. In 
any case, in the embedding space formalism we can write this integral as 

^ ^ 1 (3.1) 



^ J 2ir 2 (Pi • Q) (P 2 • Q) (P3 ■ Q) (Pa • Q) ' 
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with d 4 Q = d 4 x. We begin by following the standard route for evaluating such integrals: 
we introduce Schwinger parameters to exponentiate the propagators, obtaining 



+00 



d 4 Q 

2 - 2 Jo , , 



Y[dU e^UuP&Q. (3.2) 



After performing the Q integration (which is simply a gaussian integral), we get 



+00 



4 



i=i 



Y[dUe^=^ 2 . (3.3) 

i=i 

Up to now our manipulations have been standard, but now we do a trick, writing 

A A r „ 2 

i = dXe x T{X-Pi) (3.4) 

where M 6 is 6-dimensional Minkowski space, and the function T(X; {Pi}) is defined by: 

p+oo 4 4 

T(X;{P l })= HdtiSiX-Y^tiPi). (3.5) 

J ° i=l i=l 

The function T is known as the multivariate spline, or truncated power function. Here we 
shall simply call it the spline. Splines are best known for their applications in approximation 
theory, computer graphics, and also for the important problem of counting integer points 
in polytopes - we will comment on this latter application in the discussion. Here we see 
them making an unexpected appearance in quantum field theory! A spline is a function 
depending on a set of vectors, here the Pi, which are known as its nodes. There can be any 
number of these, but if the number of nodes is less than the dimensionality of X (in this 
case, six), the spline is not a function but a tempered distribution. To compute the spline, 
let us assume that the vectors Pi are generic, i.e. there are no linear degeneracies between 
them. Now introduce two other vectors P5, Pq such that the full set of the Pi forms a basis 
of M 6 . We can then write the spline as 

r+00 6 6 

T(X;{Pi})= l[dt i 8(t 5 )8(t 6 )S(X-J2UPi). (3.6) 

^° i=i i=i 



By thinking of the set of six Pj vectors as forming the columns of a matrix M, the set of U 
integrals is easily performed: the six dimensional delta function imposes U = Wi ■ X, with 



W M _ 1 r M 31-35 pN x pNa fo 7 \ 

Wi 5! v ,l„p/ V - v/; 7 ' - P * • (3 - 7) 

Notice that the W t M are simply the rows ofM" 1 . There is one catch though: since the t 
variables are positive, this imposes constraints on Wi ■ X. The result is then 

T(Y . P , s(w 5 -x)5(w 6 -x)Ut 1 e(w l -x) 

I {X , Pi) - , (3.8J 

yMet Pij 
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with 0(x) the Heaviside Theta function, and det Pij = det M. 



This result for the spline has a very simple and beautiful geometric interpretation. 
Going back to equation (3.5) we see that the spline can only be non-zero if X can be 
written as a positive linear combination of the vectors Pf, such positive linear combinations 
form a pointed 3 polyhedral cone, and hence the spline has support only in the interior of 
this cone. The numerator in equation (3.8) is precisely the characteristic function 4 for this 
polyhedral cone. Indeed, the W% vectors define hyperplanes on which all vectors except Pi 
lie - forming the faces of the polyhedral cone. 

Let us go back to the computation of the box integral. With our result for T we see 
that it reduces to a gaussian integral over the polyhedral cone. But we can do better! Since 
the gaussian measure only cares about the norm of X, we can do this integral straightaway. 
Indeed, since the spline is homogeneous of degree —2 in \X\ this integral can be done quite 
easily, yielding a gamma function. Actually, there is an important point: since the six- 
dimensional space has necessarily Lorentzian signature, the sign of X 2 isn't necessarily 
fixed. However, since all the Pi are future pointing, any positive linear combination of 
them will form a timelike vector, and hence X 2 < inside the cone. After the integration 
over the scale of X is made we are therefore left with an integral over vectors inside the 
polyhedral cone which further satisfy X 2 = — 1 - in other words, over an AdS slice of the 
polyhedral cone. In this case, such a slice is a hyperbolic ideal 5 tetrahedron. Let us denote 
such tetrahedron by V^. Then we have just shown that 



The X integration measure is such that the denominator above only involves the Pj with 
i = 1,...4. The box integral is equal to a volume ratio; in the denominator we have 
the volume in CP 5 of the tetrahedron with vertices denoted by the Pi, volume which is 
computed with respect to the reference vector /. In the numerator we have the volume of 
the hyperbolic "shadow" cast by that solid: the hyperbolic 3-simplex Va- This is illustrated 
in figure 2. 

For a one-loop integral such as the one we have just computed, it is quite natural to 
change the normalization so that the denominator factor changes. However we shall see 
later on that in more general computations the natural objects to consider are precisely 
these volume ratios; or put differently, AdS integrals of splines. 

4 Generalized Star-Triangle relations 

While this interpretation of the box integral as a volume in AdS is not new, and was in 
fact derived fairly recently in [33], our way of obtaining it certainly is, and sheds light 
on the mechanism at work: the essential point is that one-loop integrands are Laplace 

3 It is pointed because for all Pi we have Pf — y2, hence the origin does not belong to the cone. 
4 The characteristic function xs{ x ) °f a se t & ls defined as being unity if x G S, zero otherwise. 
5 Here ideal refers to the fact that its vertices lie on the AdS boundary. 




(3.9) 



-7- 



Figure 2. A projective triangle and its spherical shadow. A similar picture holds in the hyperbolic 
case and with tetrahedra (though harder to visualize). The projective triangle lies on the P + = 1 
hyperplane. 

transforms of splines, as we shall see in section 5. This fact immediately allows for many 
generalizations to other kinds of loop integrals. 

Let us therefore consider a more general class of one-loop integrals 6 , the colored poly- 
gons in arbitrary dimension: 

'-f&ft™^ («) 

The colouring refers to the fact that each vertex Pi now has an associated weight, Aj, 
or "colouring". We also relax the constraint that Pi is null, allowing for Pf = —rnf. 
Accordingly this implies that Pi/rrii = Xi is a vector lying on an AdS subspace of M d+2 . For 
all Aj = 1 the above integral can be understood in the dual conformal picture as a one-loop 
massive polygon integral. Non-unity values of Aj arise frequently in conformal integrals. 
Here we show that this entire class is associated with colored simplices or polytopes in 
AdSd+i- 

The computation is straightforward, and is entirely analogous to the one in the previous 
section. We introduce Schwinger parameters and perform the Q integration which is simply 
a gaussian. Generically the conformality condition E = Y17=i ^-i = d can be violated. The 
result is 

1= / Dlt(Y i t i f- d e<ZW = / dXil-Xf^e^T^MX^Xi}) (4.2) 
Jo Jm D 

6 This class of integrals has also been studied quite recently in [28] for the special case rrii = and 
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where we have introduced for convenience the measure D^t = 11^=1^^* 1 , and the 
generalized version of the spline is 

/•+OO n 

T {At} (X;{X l })= Dlt 8{X -2>Xi). (4-3) 
J ° i=l 

Notice that the spline now depends on an extra set of arguments, namely the weights 
Aj associated with each node. We will omit this argument in the case where all Aj = 
1. As before, the spline has support inside the polyhedral cone with vertices -Xj, and is 
homogeneous in |X| with degree X — D. The integral over the norm \X\ is again performed 
trivially and one is left with 

I= T{E ~ h) [ dX(I-Xf- d T {Ai} (X;{X t }) (4.4) 

^ J AdS 

After the integration over the norm the integral is now over a colored polyhedron whose 
vertices lie on the horospheres with radii m; in AdS, where again the coloring refers to 
the fact that the vertices of this polyhedron are labelled by the weights Aj. But what are 
we integrating over? To see this we need to evaluate the spline explicitly and this is done 
starting from expression (4.3). The final result depends crucially on the number n of nodes 
and on the weights Aj. We will address the cases n > D in section 6. The simplest case is 
when N = D, whereupon the spline becomes 

T{a .(*,»}) _ n&.W •*)'-' ew-*) 

Just like in the previous section, the O functions enforce that the spline is supported only 
inside the hyperbolic polyhedron V with vertices Xi, or equivalently in the interior of the 
hyperplanes defined by the Wj. The case where the number of nodes is less than D is also 
simple: we just eliminate the corresponding Wi ■ X factor and replace the Heaviside Theta 
function by a Dirac delta function, which further localizes the spline. In any event, both 
these cases lead to 

T(E-h) = r dX{I . x) z-dfr {w ■ x^-i (46) 
2^detX l .X 3 J Vn l\ y 

That is, the integral is interpreted as the computation of certain "moments" of a hyperbolic 
polyhedron V n - A direct interpretation of the integral as a volume is only possible in the 
special case S = d and all Aj = 1 - nevertheless there is still an underlying geometrical 
structure. In any case, I can still be considered as an integral of certain homogeneous 
factors over a hyperbolic polyhedron. Using the alternative form 



y/det Xi ■ X 



3 



we see that something quite interesting occurs if all dimensions Aj are integers. Indeed, in 
this case the presence of the gaussian measure allows us to integrate by parts. This gets 
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rid of the annoying X factors and at the same time the derivatives hit the functions, 
further collapsing the integral onto lower point simplices. This allows for a geometrical 
interpretation of these integrals as sums of volumes with well determined coefficients. We 
will see some concrete examples of this in section 6. 

One interesting outcome of these results, is that we see that even non-conformal in- 
tegrals (with arbitrary Aj) can suitably interpreted as integrals in AdS. But there is a 
different way of seeing it. With the parameterization, 

X = — (l,xl + x 2 ,x») (4.8) 

Xq 

we have X I = 1/xq and the AdS measure is dX = dx^ . Therefore the factors 
(X ■ I)^~ d can perhaps be better thought of as deforming the AdS measure into dx^ -4?+r> 
thereby breaking conformality. 

Let us mention another possible generalization of our results. The appearance of the 
spline can be seen to be entirely due to the factorized form of the integrand. Such factorized 
forms appear not only in flat-space integrals, but also in the computation of AdS Witten 
diagrams [26], and in particular of contact interactions, which take the form: 

The only difference from the flat-space case is that now the result is automatically con- 
formal due to the different integration measure, and hence the awkward factors of I ■ X 
never appear. We come then to the beautiful result that there is a generalized "star- 
triangle" relation in AdS space: star-like contact interactions get glued up into integrals 
over polyhedra in AdS, exchanging vertices and hyperplanes: 

HY rr r(A.) „ f . Y fr e(wyx) 

In section 7 we will see that this theme persists even beyond the one- loop/contact interac- 
tion level. 



5 Prom rational functions to geometry and back again 

In the previous sections we have seen how the factorized form of the integrand is directly 
connected to its representation in terms of the spline. In fact, we can make a stronger 
statement than this, namely that the one-loop integrand is the Laplace transform of the 
spline. Indeed, we have 

n#|^ = r°°D%te^^= I dXe^ x r {Ai} {X;{Pi}). (5.1) 
j-Jl 'Vj ' Jo JM D 

This simple fact has important consequences, since it allows us to translate identities 
between rational functions, at the level of the integrand, into geometrical identities in the 
spline representation. 
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Let us start with a trivial example. For concreteness we fix d = n = 4 and set all 
A,; = 1. Then we have the trivial identity 



Ut=iPi Q Till ^ Q 

Taking the Laplace transform this translates into 

P 5 Q 



(5.2) 



/ dXe Q - x P 5 • (-d x ) T {A} (X, {P^) 



= - j dXe- (ft • dx) ( Ilti^;Wa-^ ( 5 ,) 
Jn D \ V det Pi ■ Pj J 

where again we have introduced an extra vector Pq. The intepretation of this equation is 
straightforward: four nodes correspond to a tetrahedron; but a tetrahedron can also be 
thought of as one of the faces of a 4-simplex, which is part of its boundary. It is this fact 
that is translated mathematically by the right-hand-side of the above equation. Indeed, 
the derivative wrt X localizes the integral on the boundary of the 4-simplex defined by the 
product of theta functions and delta function. This together with the fact that P5 • Wi = 5^ 
turns the integral into one on the tetrahedron. 

If we now focus on the full integral and not just on the integrand, we have 

l q r dxe^p 5 . ( -a,)nkg(^g^) = 

J Jrd y/det Pi ■ Pj 

J 7=1 VdetPi-P? M 



where in the last step we have integrated by parts and done the integral over the scale |X|. 
In this last form, the integral is over a 4-simplex in AdS§. And yet, we know that the 
very same integral is given by the volume of the tetrahedron with vertices given by the 
first four Pi. In this way we have derived a non-trivial relation between a moment of a 
4-simplex and the volume of one of its faces. 

Now let us consider the slightly different integral, 

We will impose special kinematics: we demand that Pi ■ P^+i = 0, that is, the coordinates 
Xi and Xi+i are null separated. Further for Y we shall take the special vector 

Y = -P24P35P1 + ^14^35^2 + ^14^25^3 " ^13^25^4 + ^13^24^5 (5.6) 

It is clear that with these special choices, the five-point integral above collapses into a sum 
of five appropriately normalized box integrals. The normalization has been chosen such 
that each of these box integrals yields the volume of a tetrahedron (i.e. the det Pij is 
taken out). There is another consequence of these choices: Y ■ Pi = for all i ^ 5. Thus 
the integral above is closely reminiscent of the chiral pentagon integral [32]. 
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Its calculation is exactly the same as in our previous example and so we find 



h = Je x2 Y.(-d x )fl 

i=l 



e(Wi ■ x)5(w 6 ■ x) 

^det Pi ■ Pj 



So, this kind of integrals is associated with the computation of a certain moment of a 
4-simplex, the moment being calculated along the direction Y . 

This last form can be rewritten again. Indeed it is clear that due to our special 
kinematics, Y ■ X is only non-zero if X contains a piece proportional to P5. Such a piece 
can be captured by taking its internal product with W5, and so we also have 

15 = J*p P / dX {W5 • X) (5 - 8) 
ydet Pi ■ Pj Jv 



This result holds for any Y so far as Y ■ Pi = for all i < 5. In particular this shows that 
the choice of Y in [32] and ours lead to completely equivalent results: their Y differs from 
ours by an additional piece ~ ePi-P^Ps-Ri-Ps; the coefficient of which is chosen such that 
Y 2 = 0. Since this additional piece vanishes contracted into P5, we see that both choices 
are equivalent. Incidentally this shows that in this simple 5-point example the chirality of 
the chiral pentagon is actually completely irrelevant. 

The chiral pentagon integral, due to the presence of its special numerator factor, actu- 
ally reduces to a sum of boxes, and this can be seen directly in the spline picture. Previously 
it had been found that the n-pt 1-loop MHV amplitudes correspond geometrically to the 
volume of the boundaries of closed polytopes. The 5 pt case is an especially simple case of 
this, where this closed polytope is a 4-simplex; 4-simplices are associated with integrands 
with 5 propagators, and the presence of a numerator naturally gives the boundary of this 
4-simplex. In order for this to happen, we would need a numerator of the form Y' ■ Q, with 

y = P24P35P1 + PUP35P2 + P14P25P3 + P13P25P4 + P13P24P5 (5.9) 

which we see differs from Y by the presence of two single-mass box integrals. This of course 
reproduces the MHV amplitude representation of [32] as it should. 

Similar arguments could be made for 6-pt amplitudes in terms of a "chiral hexagon" , 
and we suspect that a slightly more complicated reasoning may well be able to shed light 
on the simple integral representations of [32] for higher point functions. 



6 Simplicial decompositions 
6.1 Method 1: Partial Fractions 

When the number of nodes is larger than D, interesting combinatorial structures appear; 
the polyhedron becomes divided into simplicial cells across which the spline has discontin- 
uous derivatives. These phenomena have been systematically studied in the case Aj = 1, 
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less so otherwise. In this section we focus on this simple case, although some of our results 
can be easily generalized. 

Recall that the definition of the /^-dimensional spline is 

r+oo 



These constraints, together with the positivity conditions ti > define a convex polytope 
in t-space. In fact, any convex polytope can be defined by such equations. It is usually 
called a variable polytope, since its shape depends on the vector X. We emphasize that 
this polytope does not have to be a simplex. Also that this polytope is completely different 
from the hyperbolic simplices we have seen before - the latter lived in "target-space" , while 
the former lives in t-space. 

To be more precise, the polytope can only exist if n, the number of nodes, is larger 
than D, otherwise the spline is a tempered distribution. For N = D the polytope reduces 
to a point, for N = D + 1 it is generically a line, and so on. In any case, it is clear that the 
spline simply computes the N — D dimensional volume of this polytope, up to a prefactor. 
In particular, for N = D, the volume is equal to one as long as X lies in the polyhedral 
cone defined by the nodes P, and we recover our previous formula. For N > D we expect 
that the volume should scale like IX^ - - , and so the spline must be given by a polynomial 
homogeneous in X of that degree. It is clear however that there should be more interesting 
structure regarding the dependence on X. Indeed, at special lower codimension surfaces in 
X-space the t-polytope can degenerate onto lower dimensional ones. Across these surfaces 
the homogeneous polynomials have to be continuous, but their derivatives do not. In target 
space, we have a polyhedral cone defined by the nodes P, and this cone is generically not 
simplicial for N > D. However, it can be decomposed as a sum of simplicial cones, which 
define a cell structure on the original cone. As X moves from one cell to the next, the 
derivatives of the polynomials can jump. 

To analyse this structure, and determine precisely the form of the spline, we can use 
the Laplace transform, or in other words, go to the integrand form that we know and love. 
There we can use elementary linear algebra to simplify the integrand and transform back 
to figure out what is happening geometrically. We begin with the scalar box integrand, 
which is the Laplace transform of a four-spline. But now, we imagine choosing kinematics 
such that the nodes Pi lie on a 3 dimensional sub-space. What this means is that we are 
picking all the Xi coordinates to lie on a one-dimensional subspace. To see that something 
must occur, recall that the scalar box integral is equal to the volume of a tetrahedron in 
hyperbolic space. If we demand that our kinematics lie on a one-dimensional subspace, our 
tetrahedron gets squashed, because it cannot live on an AdS2- This is shown geometrically 
in figure 3. Accordingly, the box integral naively vanishes. However, we are forgetting that 





We can think of the <5-function as imposing a set of linear equations in t-space, 
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Figure 3. A 4 vertex cone and its decomposition into simplicial 3-vertex ones. One needs three of 
these to obtain the same cell structure. 



actually the unnormalized box integral is given by a volume ratio: 

A — A vol Va 



(6.3) 



The denominator must also vanish for Id kinematics, since it is one of the Gram determi- 
nants, and hence we may hope for a finite final answer. This is indeed the case, as can be 
seen from the explicit expression for the box integral: 



Z. 2U 2 (-up) + 2U 2 (-vp) + log (£) i g ( ig) + log(pn) log(H + * 



with 



If we now set u 



u 



u — v) 



4uv. 



3.4) 



(6.5) 



z , v 



(1 — z) 2 + e and expand for e < 1, we find that A ~ e, and so 
the numerator of the box integral must be vanishing as well. This is of course nothing but 
the hyperbolic volume going to zero. The resulting expression is 



logfl 



log 



(-) 



log(*)-log(i) + lo § 



log 



(6.6) 



2z(l-z) 2z{\-z) 1-z 2z(l-z) 

Depending on the kinematic range - whether z is negative, between and 1 or larger than 
one - the expression can simplify. Now let us try to interpret what's going on in the spline 
picture. We start off with the box integrand 

' (6-7) 



(Pi ■ Q)(P2 ■ Q)(Ps ■ Q)(Pa ■ Q) 
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Since we are in Id kinematics, there can be at most three independent P vectors (since 
they live in the d + 2 dimensional embedding space) . Define the vectors orthogonal to the 
Pi, i.e. W} 3 -Pi = l, W 2 \ • P 2 , 3 = 0. Then we have 

A 1 W 23 • P 4 Wl-P, 

f = \ ^ Q (P2 ■ Q)(Ps • Q){Pi ■ Q) 2 (Pi ■ Q)(P* • Q)(Pa ■ Q) 2 • • Q)(Pa • Q? 

Hence the box integral has decomposed into a sum of triangle integrals, which can be easily 
done via the star triangle relation. Geometrically, the 4-vertex cone is being decomposed 
as a sum of three simplicial cones; algebraically the 4-spline is written as a sum of three 
3-splines. Consider for instance the first term: 



• Pj 

(P 2 • Q)(P3 ■ Q){Pa ■ Q) 2 " 2: ' 



23 4 -wk-a 



1 



(P 2 • Q)(P3 ■ Q){Pa ■ Q). 

= J dXe-Q x (Wl 3 -X)Ti,i,i(X,{P 2 ,P s ,P 4 }) (6.9) 

So we have a polynomial, • X living on a triangular cell. The full 4-spline can hence 
be written as the sum of three pieces, corresponding to three-overlapping simplicial cones, 
which put together make up the polyhedral cone with vertices Pi, P 2 > P3, Pa- On each piece 
there lives a different linear polynomial in X. 

This basic example can be easily generalized to a higher number of nodes or higher 
dimensionality and the conclusions are similar. So much for the computation of the spline, 
and hence for the integrand. But what does this mean for the integral? 

If we do the Q integral something very interesting happens. The e~®' x becomes e x2 , 
and hence the polynomial can be written as a derivative of that factor 

J dXe x2 (W 2 1 3 -X)T 1 , 1 , 1 (X,{P 2 ,P 3 ,P 4 }) = J dXe x2 (^ 2 V^)ri,i,i(X,{P 2 ,P 3 ,P4}) 

= - f dXe x V 2 Vd3?)7!,i,i(X,{P 2 ,P 3 ,P4}). (6.10) 

Acting on the spline we pick up the boundary of the region we are working in, so that the 
integral further localizes on the lines that make up the triangular cell! Since all the weights 
are one, we will end up getting lengths of lines in AdS? which are however divergent. Such a 
divergence can also be seen directly from the point of view of the triangle. The dimensions 
2, 1, 1 lead to a gamma function with zero argument, and a Py with zero exponent. To 
regulate we can set Pf = e and send e to zero at the end of the calculation. The length 
of a line in AdS goes like a log so we're clearly on the right track to recover (6.6). More 
concretely, denoting the length of the line between points i and j by L(i,j) we have 



log 



p. 



p p 

L(i,j) = -j=^^= (6.11) 
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and 



Id 



with e.g. 



= (Wh ■ Wi 4 ) L(3, 4) + (W} 3 ■ Wl) L(2, 4) + (W* • L(2, 3) 

+ (Wl ■ Wk) L(3, 4) + (Wl ■ Wl) L(l, 4) + (W* • W^) L(l, 3) 

+ • W24) ^(2, 4) + (Wf 2 • W 4 \) L(l, 4) + (W? 2 • W? 2 ) L(l, 2) (6.12) 



xAB pC pD p n 

A straightforward calculation now leads precisely to result (6.6). 

Although our calculations were performed for the scalar box, it is clear that the method 
is generically applicable to higher point integrals. By applying the partial fraction decom- 
position we can successively reduce any product of propagators to products involving only 
D of them. Since for such products we already know the spline representation, we are done. 
Geometrically the polyhedron gets decomposed into simplices, whose overlap reproduce the 
"big-cell" [35] decomposition of the original polyhedron. At least this is all true for the 
case where all weights equal one. To the best of our knowledge the generic case has not 
been addressed in the literature, and it indeed it seems hard to determine the spline in 
those cases, since the partial fraction certainly fails if the weights are not integers. Next 
we describe another method for determining the spline decomposition, which is also only 
valid in the unity weights case. 

6.2 Method 2: Residue computations 

We have already seen the method of partial fractions in action in the decomposition of the 
box into lines. As it turns out there is a more efficient way of doing these computations in 
terms of residues. 

Let us go back to figure 3, showing our squashed tetrahedron. Although we see from the 
picture that there are several possible simplicial cones, there is a minimal (non-unique) set 
which paves the entire cone. In this case it is easy to check that there are only three triangles 
required. Algebraically, we have a set of vectors which are not all linearly independent, 
and from these we want to look at the set of possible basis. To each basis corresponds a 
single simplicial cone. The question is then: what is a minimal set of basis one must use 
to pave the full cone? The answer is the set of unbroken basis. This is a combinatorial 
concept which can be understood as follows. Consider an ordered set of vectors Pi, and 
choose a certain basis labelled by integers i\ < %i < . . . < ip. We say that the basis is 
broken if there is a vector Pj such that j < ir> and Pj is linearly dependent on the vectors 
ifc with k > j. In other words, if this vector can be written as a linear combination of the 
basis elements to its right on the list. With this definition, it is clear that any unbroken 
basis must include the very first vector. 

In practice, if we assume generic kinematics, finding the set of unbroken basis is very 
simple. One fixes the first vector which has to be Pi, and then chooses all possible D — 1 
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subsets of the remaining n — 1 vectors. So for instance, in the four point case given 
above, a set of unbroken basis is (1, 2, 3), (1, 2, 4), (1, 3, 4), and it is easily checked that 
the corresponding simplicial cones cover the full four-point polyhedral cone. In general we 
therefore expect there to be (n— 1)1/ (D— l!)(n — D) \ unbroken basis for generic kinematics. 
So, for a 6-node spline in D = 4 (in other words, 2d kinematics since the embedding space is 
D = d+2) we expect a paving of the hexahedral cone by 10 tetrahedral cones, corresponding 
to the unbroken bases (1, 2, 3, 4), (1, 2, 3, 5), ... , (1,4, 5, 6). 

Going back to the spline, we can write it as the sum of several pieces, each correspond- 
ing to one of the elementary simplicial cones. On each of these the spline will be given by 
a polynomial, and it is therefore our next task to determine the polynomial associated to 
each unbroken basis. To determine the polynomial we compute a residue, and this is done 
as follows: if the elements of the unbroken basis are Pi k we do the map 

Pi x ■ Q ->• ui, ,P i2 -Q^ uiu 2 , etc. (6-14) 

We now take the Laplace transform of the spline, which is L = Yl™ = i(Pi ■ Q)^ 1 , do this 
map and multiply by the form 

e XQ dP h ■ Q A dP i2 -Q A ... = e x ' Q u^ l u 2 D ~ 2 . . . un^dux A du 2 A . . . A du D (6.15) 

In the above we can always decompose X as a sum of Pi so the exponential is also a function 
of the Ui. We think of these Ui as complex variables and compute the multidimensional 
residue. This is quite interesting, since we'll get an expression of the form 

i 

with f(ui) an analytic function. In terms of the original variables notice the measure is 

/d.o g (P 1 ^)d 1 „ g (^) dl o g (^|)... dl „ g (^) (M7) 

Let us consider an example, for instance when N = D + l, and let us consider the no- broken 
basis b formed by the first D elements. Then we get that 

Le X ' Q u\~ 1 u s 2 ~ 2 . . . u s ^idui A du 2 A ... Au s = 
dui A du 2 A ... A m s 1 + W[ b) ■ X Ul + W { 2 b) -Xu x u 2 + ... 

Ul-..U S Ul ( W (P) . p n + U2W (b) . p n + . . .) 

Here are the vectors satisfying ■ P^ = <5jj. The residue is simply 

wj b) • X 
wi b) -p n 

Now let us consider the slightly more complicated example n = D + 2. For each unbroken 
b let us denote P^^P^ the two elements which are not in it. Then it is easy to check 
that the polynomial corresponding to a given basis b is simply 

{W[ b) ■ X) 2 

{w[ b) ■ P[ b) ){w[ b) ■ P?) 



(6.18) 



(6.19) 



(6.20) 
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This means that the full spline can be written in this case as 



T(X;{P l }) = Y^ 



(6.21) 



{6} 



where \(b) ls the characteristic function of the simplicial cone corresponding to the basis 
b, namely 



D 




i=l 



and det b T b is the determinant of the matrix of the Pij, where the i, j run over the elements 
of the basis b. 

Again, we have illustrated the residue method with these simple examples, but more 
general cases are also easily determined. In the non-generic case where there are extra linear 
deneracies the algorithm is unmodified; one must simply be more careful in constructing 
the set of unbroken basis. 

6.3 Doing the integrals 

Now we come to the interesting part. We are not interested in the spline, but rather in its 
integral over a gaussian measure, 



Let us consider the case n = D + 2. Then the integral above is a sum of several pieces, 
each of which takes the form 



But now we see something very interesting: because of the presence of the gaussian measure, 
we can write the polynomials as derivatives, and integrate by parts! There are then two 
kinds of terms that appear. One of them is simply 



which is, up to an overall factor, the volume of the D — 1-dimensional simplex in AdS 
with vertices given by the unbroken basis b. The other kind of terms correspond to lower 
dimensional simplices, which arise by hitting the characteristic function with derivatives. 
Now, the remarkable thing is that if the all simplices are ideal (which means that all the 
Pi square to zero 7 ), these terms seem to cancel among themselves! Although we have 
no formal proof, we have checked in several examples that this is the case. One way of 

7 In terms of loop integrals, this means that all propagators are those of massless particles. 




(6.23) 




(6.24) 




(6.25) 
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understanding this might be that such terms would lead to integrals with necessarily lower 
transcendentality than those such as 6.25. 

To summarize, we conjecture the following: the d-dimensional ci-gon loop-integral can 
be written, upon using d — A dimensional kinematics as a sum of (d — 2) dimensional d — 2- 
gon loop integrals. In particular this implies that an hexagon integral in 2d kinematics 
can be written as a sum of 10 box integrals; and that an octagon integral in Ad kinematics 
can be written as a sum of 7!/5!2! = 21 hexagon integrals. Each integral is labelled by an 
unbroken basis containing d — 2 elements, and the coefficient of each such integral is 

(w? ■ P?)(w[ b) ■ P?) 1 J 

This is only one of many other possible conjectures which would hold for more complicated 
integrals with N > D + 2. Work is in progress on these matters [38]. 



7 Spline convolutions and multiloop results 
7.1 Splitting the tetrahedron 

Let us go back one last time to the four-dimensional scalar box integral. Thanks to the 
spline, it can be nicely interpreted as the 3-volume of a hyperbolic tetrahedron (up to a 
normalization factor). Now, it is clear that at least in Euclidean space, one can also see a 
3-volume as roughly the product of a length by an area; or as a triple product of lengths. 
So is there an analogous statement in hyperbolic space? We shall now see that indeed there 
is, and the answer is connected to the circumstances under which splines can be merged 
together. Indeed, let us again consider the box integrand: 

(7.1) 



(Pi • Q) (P2 ■ Q) (Ps ■ Q) (Pa ■ Q) 

The connection of this integrand with the four-noded spline is directly related to the pro- 
cedure of joining the four separate propagators into a single one by use of Feynman param- 
eters. But this is not the only choice: we could for instance join three of them, or make 
two groups of two. This leads to the appearance of several splines, which however must be 
exactly equivalent to the original one. We can see this quite directly. Recall that we have 
found 



r r r+oa 4 

/ dXT(X,{P u A t }) = dX Y[dtiS(X-Y,iiPi)- (7-2) 

^AdS J AdS JQ „•_ i 



1=1 

Now we can repeat our original trick of introducing delta functions. For instance, 

" + 00 



r f r+oo ' 

/ dXidX 2 / dX / Y[dti5{X-Xi-X 2 ) 

Jm d Ads Jo i=1 



5(X 1 -t 1 P 1 -t 2 P 2 )5(X 2 -t 3 P 3 -t 4 P 4 ). (7.3) 
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To continue, let us separate the integrals over the scale factors |-Xi|, \X 2 \ and denote these 
two variables by s\, s 2 . Then it is easy to find 



4 



= [ dXdX ± dX 2 f °° ^^sjsl 5{X - 8l X! - s 2 X 2 ) x 
Ja&s Jo s i «2 

/■+00 r+00 

/ dt 1 dt 2 S(X 1 -t 1 P 1 -t 2 P2) / dt 3 dt A 5{X 2 -t 3 P 3 -t A P A ) 
Jo Jo 

= [ dXdX 1 dX 2 T2,2(X,{X 1 ,X 2 })T(X 1 ,{P 1 ,Pt i })T{X2,{P^,PA}) (7.4) 

^AdS 

Written in this form this might seem a bit obtuse, but the interpretation is actually quite 
simple. Firstly, there are the two splines with arguments X\, X 2 . These have support 
on the hyperbolic geodesic lines from P\ to P 2 and P 3 to P4 respectively. The remaining 
spline has support on the hyperbolic geodesic connecting X\ and X 2 itself. Notice though 
that the weights of the nodes are no longer 1. This seems to be the closest one can get 
to a natural splitting of the volume into three "lengths". The use of inverted commas 
is very important here, since as we've shown before, the integral over the spline gives at 
best a volume ratio, and only in the case where all the weights are equal to one. In this 
case we see that there is no getting around the fact that the spline is a ratio; no factor 
can be multiplied to cancel the denominator of the spline with argument X, since such a 
denominator involves X\,X2 which are being integrated over. 

It is also easy to show that there is an alternative representation in terms of only two 
splines, 



I dXdX' Ti, 3 (X, {P 4 , X'}) T{X', {Pi,P 2 , P 3 }) 

./AdS 



dX 1 T(X, {P 1 ,P 2 , P 3 }) (7.5) 

AdS 1/4 " I 

Notice the appearance of l/(Pi • X), which is nothing but the scalar bulk to boundary 
propagator in AdS\ Such a factor is precisely the AdS integral of the spline with nodes 
Pa,X. Since the bulk-to-boundary propagator is a limit of the bulk-to-bulk propagator, 
one may naturally wonder what the result would be if instead of the massless P4 we 
had a massive, AdS vector. We show in the appendix that J Adg dX 7ai,a 2 (^ {-^Ij X2, }) 
is essentially the sum of two AdS bulk to bulk scalar propagators. In the case where 
Ai = A2 = 2 = d/2 one obtains the logarithmic bulk-to-bulk propagator of a field of 
dimension d/2 in AdS. The conclusion is that the original flat space integral is given 
entirely in terms of natural quantities in AdS space. 

There are several obvious generalizations of these results. The number of nodes in- 
volved and their type (massive or not), the dimensionality of spacetime, the weights and so 
on are all quite irrelevant. In this way, e.g. a hexagon integral can be written in a number 
of different ways: as the volume of a hyperbolic 5-simplex; two hyperbolic triangles glued 
by a line; a triangle with one of its tips lying inside a tetrahedron; and so forth. The rule is 
that if a given node X' of a spline T' is the argument of a spline T", and if this argument 
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X' is being integrated over, then both splines can be merged into a single one as long as 
the sum of the weights in T" is equal to the weight of the node X' in T' ■ So, for instance 

f dX'Ti,3(X,{P 4 ,X'})T(X',{P 1 ,P 2l P 3 }) = T(X,{P 1 ,P 2 ,P 3 ,P 4 }). (7.6) 

J AdS 

These rules are useful since sometimes a given representation in terms of some number of 
splines might be easier to derive than others, especially when one is considering multi-loop 
integrals. It is to these that we turn next. 

7.2 Multiloops 

For definiteness, we shall start by considering the fully-massive double-box integral, which 
in the dual conformal picture corresponds to a 3-to-3 tree-level exchange diagram. This is 
displayed in figure 1. The integral is written 



V 




= 


A 





dQidQ 2 

(2vr rf / 2 ) 2 Qi • Q 2 i\ Pi ■ Qi i\ Pi ■ Qi 



Our first question is: can such a integral be written solely in terms of AdS quantities, or 
was this rewriting a lucky accident at one-loop? As usual, introduce Schwinger parameters 

J ^ + °° ds JJa*. p df . e Qi'(E, 6i ^)+Q 2 -(E ie H*^i)+^i-Q 2 (7.8) 

There are several ways to proceed now. We could introduce splines, one for each side of the 
diagram. This would correspond to two triangles connected by a convolving kernel which 
turns out to be the integral of a 2-spline. Instead, we first do the Q\ integral, giving 

J ^ d ' 2 ) Jo i L L W 

Performing the Q 2 integral we can write the result in turns of two splines, 

f dX 1 dX 2 e x '+ x 'T(X 1 ,P L )T(X 2 ,{X 1 ,P R }). (7.10) 

JM 6 xlM« 

with Pl,r = {Pi,i £ L,R}. Homogeneity in the norms of X\ and X 2 allows us to perform 
these integrals easily, and we find 



X 



I dX 1 dX 2 T(X u P L )T(X 2 ,{X 1 ,P R }) (7.11) 

J AdS 



The interpretation is clear: we are doing an integral over a tetrahedron with vertices 
X\, P4, P5 and Pq; at the same time we are integrating over the position of vertex Xi, but 
not over all of AdS, rather inside the hyperbolic triangle with vertices Pi,P2,P3. 
Using our merger rules, we can derive the alternative representation 



fTl jAdS i=i 



3 

/ Yl dX i T ( X ^ p L)TiM^2,{Xi,X 3 })T(X 3 ,P R ). (7.12) 

J AdS „•_ 1 
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Figure 4. An AdS picture of the double box integral. 



which can be visualized geometrically in figure 4. We have split the tetrahedron into a line 
and a triangle, so that the very same integral is seen as two hyperbolic triangles connected 
by a line. More precisely, this line is precisely the two-spline mentioned in the previous 
section, and computed more carefully in the appendix, corresponding to the sum of two 
bulk-to-bulk scalar propagators in AdS; one each of dimensions 3 and 1(= d— 3). We could 
of course keep going and further split the triangle into a pair of lines. However, it seems 
impossible to join the tetrahedron and triangle into a single object. 

Of course, everything we have said up to now is easily generalizable to more complicated 
integrals. Weights can be varied, masses can be added, and so on. Once this is done, 
another kind of phenomenon appears. Consider a 2 — > 2 exchange diagram with given 
weights Aj and an internal propagator of weight 5. Then it is not hard to show that this 
diagram is given by 

r ( Al + f 2 "* ) / pd* rw^^xr) 

Ts^-s(Xi,X 2 ] X 3 ) Ta 3 ,aAP3, Pa; X 2 ). (7.13) 

This is quite interesting, because the 4-pt contact interaction admits a very similar repre- 
sentation, once we split the tetrahedron onto a product of three lines: 

r 3 

/ 11^1 7ai,A 2 (Pi ; P2;^i)7ai+A2,A3+A 4 (^1>^2;^3)'7a3,A4(- P 3,-P4;^2)- (7.14) 
i=l 

Equality can be only be achieved if Ai + A2 = 5, at which point the exchange diagram 
develops a pole from the gamma function. But so, we have learned that by looking at 
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Figure 5. A (fully massive) 4-loop integral, its position space dual and its glued up version in AdS 
space. 



the behaviour of integrals when the weights are varied, we can find new connections be- 
tween integrals with different numbers of loops. In the case in point, we see that up to a 
controllable divergent pre-factor, the one-loop integral is a special case of the two-loop one. 

Let us close this section by a few further comments. This two-loop example shows that 
the generalized star-triangle relations continue to hold beyond one-loop. Stars get glued up 
onto hyperbolic simplices, and flat-space integrals transform into AdS ones. For instance, 
the symmetric 12-pt diagram in (j) 4 theory gets beautifully glued up onto a picture of a 
hyperbolic tetrahedron whose vertices lie on ideal hyperbolic triangles, as shown in figure 
5. 

Many other different pictures are possible by splitting the tetrahedra in various ways, 
or joining some of them onto higher dimensional polytopes. Also, by appropriately vary- 
ing the dimensions of the internal propagators poles can develop which relate this integral 
onto simpler ones; such as the 12-pt star. Therefore this part of the story also general- 
izes. Finally, one may well wonder about position-space loop integrals. These match onto 
"window" -type diagrams in the original momentum space. A systematic exploration of 
these diagrams is beyond our reach, but we can consider the particular example given in 
figure 6. The position space integral takes the form 

* dQj A r(AQ T(5 12 ) T(5 13 ) T(5 23 ) 

ri 2nd/2 ri (p * • ft)Al (Qi • Q2)Sl2 {Qi • Q3)8ui (Q2 • Qs)523 ' 

The trick now is to realize that 

r(s 12 ) r(5 13 ) v(s 23 ) r AV A m) 



(Qi ■ Q 2 ) 5 - (Qi • Q 3 ) s ^ (Q 2 • Q 3 ) 5 - A dS dX 2 (Qi ■ xfi (7 - 16) 

with 5\ 2 = 5\ + 5 2 — 5 3 , and cyclic permutations thereof. With the three Qi variables so 
decoupled, we can use our usual tricks to turn the corresponding integrations into AdS 
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(a) (b) 



Figure 6. A position space loop integral (a) and its AdS counterpart (b). 

ones. In this case the result is that the triangle gets "unglued" onto a star interaction. 
Since we only included one external point per vertex what we are seeing is three AdS lines 
meeting at a common, integrated point. Had we included more external points, we would 
have obtained simplices joined together by one common vertex. This result suggests that 
perhaps other position-space loop integrals can be similarly unglued onto AdS interactions. 
However, the story will have to be necessarily more complicated since already at the 4-pt 
level there is no obvious AdS amplitude which yields the right product of four propagators. 

8 Discussion 

In this note we have introduced spline technology to quantum field theory. More precisely, 
we have shown that exact calculations can be rewritten in terms of splines. The spline seems 
to act as a bridge between the integrand and its geometrical interpretation, translating 
identities on one side to the other and vice-versa. We believe that the spline representation 
will be useful in establishing the dictionnary between integrand and geometry originated 
in the works of [29, 30]. 

This dictionnary works not only for the fashionable dual-conformal loop integrals of 
N = 4 SYM, but in fact for any planar loop-type integral, with varying weights, masses 
or no masses, numerators, and so on. All these become associated to a certain colored 
polytope in AdS. Beyond one-loop level an intriguing structure appears, whose meaning 
has so far eluded us - one obtains polytopes with support on other polytopes. However we 
believe that the "charge" interpretation of the spline, to be discussed shortly, will play an 
important role in elucidating this structure. 
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With these results we have seen that a very large class of flat-space integrals can be 
completely translated into a dual computation in AdS space. But even correlation functions 
computed directly in AdS using Witten diagrams also fall under the charms of the spline, 
and we get a beautiful correspondence between star integrals in AdS and their glued up 
versions as hyperbolic polyhedra. We should note it is not the first time that a connection 
between Feynman integrals and hyperbolic polytopes has been made. In [46] this same 
point was made for one-loop integrals. Their approach can be summarized and related to 
ours in the following way. If we begin with the one-loop expression we can write 



From this point of view it is clear the value of the integral is the same as the volume of 
a certain ellipsoid in Schwinger parameter space. Depending on the eigenvalues of the 
metric Pij this can be a spherical or hyperbolic ellipsoid. We can play the same trick with 
multiloop integrals using our spline expressions but the resulting expressions then no longer 
have a simple interpretation. 

Besides the theoretical and conceptual questions posed by the existence of the spline 
representation, we have also endeavoured to show that there are very practical and useful 
consequences of this representation. One of these regards simple algebraic decompositions 
of higher point integrals in lower dimensional kinematics, a theme which is currently un- 
der investigation [38]. Another, is the simple derivation of identities between conformal 
integrals with differing conformal weights. Indeed, augmenting the conformal weight by 
unity introduces a factor W ■ X in the spline. Since the splines are being integrated with a 
gaussian measure, we can trade this factor for a derivative, thereby establishing relations 
between integrals with differing conformal weights. For instance, there are many different 
relations between conformal blocks [47-49] with different weights, spin and so on. It is 
our belief that some of these relations can be rederived from the spline picture. But we 
do not have to stop there; any one-loop integral will satisfy similar relations which can be 
explicitly derived from the spline picture. For instance, it is straightforward to check how 
this works for the 2-^1 function, which is related to the 2-node spline, as shown in appendix 
B, and one indeed recovers the usual transformation properties of the 2-^1 upon integral 
shifts. It might be interesting to work out some of these rules for more complicated higher 
point integrals. 

Splines have an important application that, though thus far has been unmentioned, 
suggests some ideas on its meaning. This application refers to the counting of integer 
points in polytopes. For physicists, this is an important problem which we can reformulate 
as follows: in how many ways can a given charge Q be decomposed in terms of a set of 
elementary charges q{! That is, can we compute the degeneracy, defined by 
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Notice that the charges do not have to be scalars. It is clear that the above is simply 
a discrete version of the spline, and it counts the number of integer points inside the 
polytope defined by M ■ T = Q, with M = (q±, . . . q n ) the matrix whose columns are the 
charge vectors. For large enough Q, the number of such points is well approximated by 
the volume of the polytope itself, and the degeneracy becomes precisely the spline! 

r+oo n 

d {qi} (Q)^T(Q,{q l })= JI^W ( 8 - 3 ) 

•'° i=l 

In other words, the spline is the continuous limit of the degeneracy, and in fact it is 
its leading approximation for large Q. Now, how would one go about computing the 
degeneracy? Physicists have also understood how to do this: one computes a partition 
function introducing a chemical potential fi, 

n 

w = e = n (8 - 4) 

Q i=l 

and reads off the degeneracy by doing a power expansion. But now we see that in the con- 
tinuous limit, forming the partition function corresponds simply to computing the Laplace 
transform, and we know that the transform of the spline is the integrand! And indeed, if 
we expand the partition function above for small qt we obtain a product of "propagators" : 

n 

zw^n— . (as) 

With this information, what meaning can be given to the integral? We have seen that a 
one-loop integral can be written as 

1=1 dXe x2 T(X;{Pi}). (8.6) 
JM d + 1 ' 1 



In our discrete analogy this would be 

Q 



This has again the flavour of a partition function. This is clear if we think of the elementary 
charges qi as momenta, since then we can think of Q 2 as energy and we have 

I = Tre~ H (8.8) 

This is not a number, but rather a function of the elementary charges, which we can think 
of as implementing restrictions on the trace. These results are suggestive that there is a 
deeper structure behind the spline representation of the loop integrals. For now we leave 
further investigations for the future. 

To end on another speculative note, let us mention the curious similarities between the 
splines and the twistor space delta functions which are the bread and butter of iV = 4 SYM 
amplitudes calculations [50]. Perhaps the latter could be understood as a complex variable 
generalization of splines. Then the geometrical interpretation of scattering amplitudes in 
N = 4 SYM would be one step closer to being a reality. 
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A Relationship to Carlson's Dirichlet B-splines 

From the Dirichlet spline or generalized truncated power, 

r {Ai} (x,{p t })= / n-r^^-E*^) ^ 

JO H 

we can go to Carlson's Dirichlet B-spline [36], by removing the overall scale. More precisely 
the T-spline is supported on a polyhedral cone, so by "gauge fixing" the scale, we will get 
an object supported on a polyhedron. This can be done by introducing a partition of unity 

Jdv5(v-Y,ti) (A.2) 
followed by a rescaling of the U, ti — > vti. The one loop integral becomes 

/ = J dXe x2 T {Ai} (X;{P l }) = J g£ {Aj} (X;{P i} ) (A.3) 

with 

B {Ai} (X;{Pi})= [ fl^^a-E^^-E^) ^ 

J i=l % 

Intuitively, going from T to B corresponds to taking a flat slice of the polyhedral cone. 
Indeed, B has support in the convex hull of the vertices P,, or in other words, the polyhedron 
with vertices Pj. Carlson introduces the concept of Dirichlet average F of a function /, 
which can be formulated as the convolution of a function f(X) with the spline. More 
precisely, we have: 

F(A i ,P i ) = J d D Xf(X)B {A} (X,{P t }) (A.5) 

The one-loop integral is then the Dirichlet average of a power function, which is one of the 
better studied examples. Carlson defines the resulting function as an P-hypergeometric 
function. Some elementary properties of the P-spline which pass straightforwardly onto 
the T-spline is that when knots coalesce, they can be replaced by a single one, with a 
weight corresponding to the sum of the weights of the previous ones, e.g. 

P({Ai, A 2 , . . . , A n }, {P l5 P l5 . . . , P n }) = F({Ai + A 2 , . . . , A n }, {P 1; . . . , P n }) (A.6) 

n n n—l n— 1 

Also, by analytic continuation, the requirement that the Aj have positive real part can be 
relaxed. In this case a knot with zero weight may be omitted. 
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B A two-loop computation 

We start off with the double box integral, or exchange diagram in position space 

[ d d Q 1( i d Q 2 T(5) j j T(A t ) yj r(Aj) 
^m-J (27rd/2)2 { _ Qi .q 2) 11 (_ P . . g l} 11 ( _ P . . g 2 ) ^ 

Introducing Schwinger parameters and splines we can write 
I n , m = J MD dX 1 dX 2 Ts L (X 1 )Ts R (X 2 ) J ^ i g^jr +0 °^/ e Q 1 -x 1+ Q,x 2+s Q 1 .Q eB2) 

where 7~s(X) denote the splines for the simplices S, which are described not only by their 
vertices Pi but also by a coloring of these vertices given by the weights Aj. To proceed we 
first notice that the overall scales |Xi| and \X 2 \ can be decoupled from the splines; since the 
splines have support on the cone formed by the null vectors Pi, the integrals have support 
only where Xf < 0, so this is well defined. Indeed we have 

T Sl {Xi) = IXtf^TsAXi) = |Xi r^TsjXx) (B.3) 

with Xf = —1, and where we used the conformality condition Y>l + 5 = d = D — 2. Defining 
the integral over the coloured hyperbolic simplex by 

[ dX= [ dXS(X 2 + l)Ts(X) (B.4) 
Js Jm d 

Then we can write 



I n , m = f dX t f dX 2 G s {X l ,X 2 ) (B.5) 
with the convolving kernel 

G S (X 1 ,X 2 )= [ +0 ° ^p^(v lV2 r s / [ ^g| e ^Q 1 ^ + 2, 2 Q 1 .x 2+2s Q,Q 2(B6) 
J Q vi dv 2 s J (27r d / 2 y 
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Notice that at this point we have reduced the computation of I n ,m to a double integral 
over AdS space. It is our next job to figure out what this convolving kernel is. Performing 
first the Q\ and Q2 integrals we find 

G S (X 1 ,X 2 ) = [ + °° ^^^(v^f-s s s e - ( i +s >?-„! +2s ^ 2 x,x 2 (B7) 

Jq Vl dV2 s 

Doing rescalings s — > s/v\ the v\ integral is simply performed. Renaming V2 — > s and 
defining the chordal distance u = (X\ — X2) 2 we find 

G 5 (X 1 ,X 2 ) = [ +0 ° — ^_ s S s d-S e -(s + sr-ssu (B8) 

2 Jo s s 

This is remarkably close to the expression for the bulk-to-bulk propagator of a scalar of 
dimension 5 in AdSd+i- Notice that this expression follows from 

(d-26\ f dQ TjS) Tjd - 5) 
G S (X 1 ,X 2 ) - r ^— j J - Q)6 - Q)d _ s (B.9) 



Continuing, 



/'dr., _,. ds ds s s_ Csd „s-c e -(s + s)i~ssu 



= r 



s s 

d-25\ I dc r(c)r(d/2 - c)r(<5 - c)F(d - S - c) 



2 J J 2m T(d - 2c) 

¥ (d-25\ I dc r(c)r(<5 - c)T(d - S - c) fu\-c 



2 d ~ l V 2 / 7 2vri Y{\ + I - c) 

r(^)r(«y)r(d-«y) p / id u 



Now we can use a hypergeometric identity: 

2 F 1 (a,b,c,z) = (^—^4) {-z)~ a 2 Fi(a,l + o-c,l + a-6,-) 

Vr(6)r(c-a); « 

which holds whenever 6 / a + n for all n € Z, and z ^ (0, 1). To us this imposes that 
5 7^ d/2 + n. When this is true, we can use this identity to find 

=, /d-2<5\ f / r(d-2^rffl \ / r(2^-d)r(d-<5) \ 1 

G S (u) = ^zi r 1^^— J 1 1 r ^ i+d-2^ - 1 Gs w + I — r( 1+2 2 g ~ d ) — ) ( ^ J 

= ^ r (^) { 2— r (^) r (*) + 2— r (^) r (d - *) G(W („) 

where 

G A (u) = 

w- 

is, up to a normalization factor, the bulk-to-bulk propagator for a field of dimension A. 
This shows that the double box integral has a nice geometric interpretation: it corresponds 
to an integral over two simplices glued by the exchange of scalars propagating in AdS. 



G A (n) = ^ 2 F 1 [A, 1 + 2 f d ,l + 2A-d,--) (B.12) 
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